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Abstract: We give arguments for the existence of exact travelling- wave solutions, of 
old (in particular solitonic) and new type, of a perturbed sine-Gordon equation on the 
real line or on the circle, and classify them. The perturbation of the equation consists 
of a constant forcing term and a linear dissipative term. Such solutions are allowed 
exactly by the energy balance of these terms, and can be observed experimentally e.g. 
in the Josephson effect in the theory of superconductors, which is one of the physical 
phenomena described by the equation. 
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1 Introduction 

The "perturbed" sine-Gordon equation 

Vtt ~ fxx + sin ip + aipt +7 = 0, x £R, (1) 

(a > 0, 7 £ R are constants) has been used to describe with a good approxima- 
tion a number of interesting physical phenomena, notably Josephson effect in 
the theory of superconductors f5], or more recently also the propagation of local- 
ized magnetohydrodynamic modes in plasma physics 20 . The last two terms 
are respectively a dissipative and a forcing one; the (unperturbed) sine-Gordon 
equation is obtained by setting them equal to zero. 

In the Josephson effect ip(x,t) is the phase difference of the macroscopic 
quantum wave-functions describing the Bose-Einstein condensates of Cooper 
pairs in two superconductors separated by a very thin, narrow and long dielec- 
tric (a socalled "Josephson junction"). The 7 term is the (external) socalled 
"bias current", providing energy to the system, whereas the dissipative term 
aipt is due to Joule effect of the residual normal current across the junction, 
due to electrons not paired in Cooper pairs. Additional terms can be added to 
describe additional features, e.g. a term like aip x would approximately describe 
|15j the Josephson effect in a junction with a linearly varying (albeit very small) 
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breadth, whereas a dissipative term like —£<f X xt would take into account (see 
e.g. [2j) the effect of the residual normal current along the junction. We plan to 
address the latter, third order equation elsewhere, exploiting our results [3 E] 
about equations involving the differential operator —ed^dt + d% — d%. 

Among the solutions of the sine-Gordon equation the solitonic ones are 
particularly important, in that they describe stable waves propagating along 
the x-line. There are strong experimental (see e.g. [2] for the Josephson ef- 
fect), numerical jSJ EI and analytical indications that solutions of this 
kind are deformed, but survive in the perturbed case; nevertheless, up to our 
knowledge there is no rigorous proof of this. The analytical indications are ob- 
tained within the by now standard perturbative method \W\ I12j based 
on modulations of the unperturbed (multi)soliton solutions with slowly vary- 
ing parameters (tipically velocity, space/time phases, etc.) and small radiation 
components. (This is inspired by the Inverse Scattering Method). The Ansatz 
for the approximate one-(anti)soliton solution reads 

<p(x,t) = g (x-x (t)-v(t)t\ +e<pi(x,t)+... 

where (fo(x,t) := go(x — vt) is one of the unperturbed (anti)soliton solutions 
given below in (|14|). whereas the slowly varying xo(t),v(t) and the perturbative 
"radiative" corrections e(pi(x,t) + ... have to be computed perturbatively in 
terms of the perturbation ef of the sine-Gordon equation (in the present case 
one may choose e = 7 and ef = +aipt + j). One finds in particular approximate 
solutions with constant velocity 

v(t) = Voo :=±[l + (4a/7r 7 ) 2 ]-5 (2) 

which are characterized by a power balance between the dissipative term atpt 
and the external force term 7. They are interpreted as approximating expected 
exact (anti)soliton solution. The experimentally observed velocity is consistent 
with the value Voo within present experimental errors. 

The purpose of this work is to give non-perturbative arguments for the 
existence of exact travelling-wave (in particular solitonic) solutions of the above 
equation on the real line or on the circle, and a preliminary classification of 
them. Our approach is less ambitious, in that it is based on pushing forward 
the study of the ordinary differential equation which is obtained by replacing 
in the equation (pQ) the standard travelling-wave Ansatz 

<p(x,t) = g(€) = g(x - vt) (3) 

(here and in the sequel £ := x — vt), and therefore cannot be applied to mul- 
tisolitonic solutions. The ordinary differential equation is the same as the one 
describing the motion along a line of a particle subject to a "washboard" po- 
tential and immersed in a linearly viscous fluid, and therefore the problem is 
essentially reduced to studying this simpler mechanical analog where £ plays 
the role of 'time'. After reviewing (sectional) travelling-wave solutions of the 
sine-Gordon equation in section |3] we classify the possible solutions of the lat- 
ter, identifying those having bounded energy density at infinity and possibly 
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yielding stable solutions for (^Q); the main results are collected in Proposition 
^ The actual existence of such solutions is strongly suggested by our physical 
expectations on the mentioned particle mechanical analog and will be proved 
elsewhere together with some general properties of these solutions (there we 
will also provide an alternative perturbative scheme for their determination). 
Among the solutions there are: those yielding solitonic solutions for (fTJl. which 
are characterized by their going to two neighbouring local maxima of the poten- 
tial energy of the particle as £ — > ±00; those yielding "array of (anti)solitons" 
solutions; those yielding "half-array of (anti)solitons" solutions. The occurrence 
of the latter is a new phenomenon, with no counterpart in the pure sine-Gordon 
case. Also, contrary to the unperturbed case, the propagation velocity v of the 
soliton turns out to be not a free parameter, but a function of 0,7, which 
coincides [H], at lowest order in 7, with ©. 



1.1 Preliminary considerations: 

Given a solution one can find a two-parameter family with infinitely many 
others by space or time translations; therefore, for each family it suffices to give 
just one representative element. 
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Figure 1: Mechanical model for the sine-Gordon equation, (a) Spring, (b) 
dolder, (c) brass, (d) tap and thread, (e) wire, (f) nail, (g) and (h) ball bearings, 
(i) base (After A. C. Scott ^5], courtesy of A. Barone, see |2) 



A system obeying the sine-Gordon equation can be modelled by the 
discretized mechanical analog in fig. ^ namely a chain of heavy pendula con- 
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strained to rotate around a horizontal axis and coupled to each other through 
a spring applying an elastic torque; one can model also the dissipative term 
—oupt °f © by immersing the pendula in a linearly viscous fluid, and the forc- 
ing term 7 by assuming that there is a uniform friction between the pendula 
and the horizontal axis, and that the latter rotates with constant velocity. 

The constant solutions of Q are ip(x,t) = — sin -1 7 + 2irk and tp(x,t) = 
sin -1 7+(2A;+l)7r. The former are stable, the latter unstable. To see this one 
just needs to note that they yield respectively local minima and maxima of the 
energy density 



This is visualized in the mechanical analog in fig. ^ respectively by config- 
urations with all pendula hanging down or standing up. We choose the free 
constant K £ M so that it gives a zero energy density at one (particular) sta- 
ble constant solution, tp(x,t) = — sin -1 7: then K = yl — 7 2 + 7 sin -1 7. In 
general, as a consequence of h fulfills the equation 



where we have introduced the energy current density j := <p x <Pt- If a = this 
is a continuity equation. The negative sign at the rhs shows the dissipative 
character of the time derivative term in if a > 0. 

Our working definition of a (multi)solitonic solution ip is: A) It is a stable 
solution which significantly differs from some minima of the energy density h 
only in some localized regions; this means that mod. 2n it must be 

lim (f(x, t) = — sin _1 7, lim f(x, t) = — sin _1 7 + 2nir (6) 

X'^ — 00 x— » + oo 

with n G 7L. Moreover, as usual we require that: B) In the far past cp is approx- 
imately a superposition of single solitons, antisolitons (and possibly breathers), 
which in the far future emerge from collisions again with the same shape and 
velocities with which they entered [however, in this work we shall only deal with 
single (anti) solitons]. As we shall recall below, one-soliton and one-antisoliton 
solutions are characterized by n = 1, — 1 respectively [whereas the static stable 
and constant solution tp{x,t) = smT 1 ^ (mod 2ir) correspond to n = 0]. In the 
mentioned mechanical analog the one-(anti)soliton solution describes a localized 
twisting of the pendula chain by 2ir (anti)clockwise, as depicted in figure|31 (a), 
moving with constant velocity. The above condition yields an energy density h 
(rapidly) going to 0, 2n7T7 respectively as x — > —00, 00; only if 7 = 0, i.e. if the 
values of the potential energy at all the minima coincide, the potential energy 
h vanishes at both x — > —00,00, and one recovers the standard definition of 
solitons. 

Although n 7^ makes the total Hamiltonian 



h '■= Y + -y + 79? - cos ip + K. 



(4) 



d t h - d x j = -aipf. 



(5) 




(7) 



—00 
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divergent, it gives a well-defined, non-positive time-derivative 

oo 

H = - j cc(f$< 0, 

— oo 

as a result of integration of (J5J). The effect of 7 7^ is to make the values 
of the energy potential at any two minima different; this leaves room for an 
indefinite compensation of the dissipative power loss by a falling down in the 
total potential energy, and so may account for solutions not being damped to 
constants as t — > 00. 

Without loss of generality we can assume 7 > 0. If originally this is not the 
case, one just needs to replace p> — > —ip. If 7 > 1 no solutions cp having finite 
limits and vanishing derivatives for x — > ±00 can exist, in particular no static 
solutions. If 7 = 1 the only static solution 99 having for x — ► ±00 finite limits 
and vanishing derivatives is ip = —ir/2 (mod 2tt), which however is unstable. 
In the sequel we shall assume < 7 < 1. 



2 The sine-Gordon equation 

Before going on, let us recall ^ |3] what happens in the case 7 = = a. The 
sine- Gordon equation is 

<Ptt ~ <Pxx + sin p = 0, x£t. (8) 

The associated Hamiltonan (a conserved quantity) is 

H := J f^ + ^-cosp + lW (9) 

—00 

Let £ := x — vt. One looks for travelling- wave solutions of ((HJ), namely for 
solutions with a 'fixed-profile' Ansatz ©. As the time t flows the profile of g 
will move to the right/left according to a positive/negative sign of v. Replacing 
the Ansatz in the sine-Gordon equation one finds for g the equation (v 2 — X)g" + 
sing = 0. If v 2 = 1 the equation admits only the (stable or unstable) static 
solutions ip(x,t) = g(£) = 0, tt (mod 2tt). Otherwise it can be rewritten in the 
form 

g" + smg = 0, (10) 

where 

where rj,a = ±1 are arbitrary signs. This is the equation of motion w.r.t. the 
'time' £ of a pendulum {g being the deviation angle from the stable equilibrium 
position), or equivalently of a particle in a sinusoidal "energy potential" U = 
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Figure 2: The potential energy U 



— cos g. Multiplying the equation by g' one finds the "mechanical energy of the 
pendulum" integral of motion e 



d_ 



g' 2 

cos q 

2 y 



q' 2 

e := + U = const. (12) 



By definition e > — 1. The above implies 

±dg 1 = d£. 

y / 2(e+cosg) 

The latter equation can be immediately integrated out, giving £ = £(g). The 
final step is the inversion of this function in disjoint intervals, which will give 
g = <?(£), and patching the intervals. 

According to the choice of e, v one obtains different kinds of solutions. Plot- 
ting the potential energy (Fig. [2J helps us get an immediate qualitative under- 
standing of them. 

1. If e = eo := —1 then necessarily g(£) = 0(mod27r). This corresponds to 
the constant solutions 



ip s (x, t)=0 (mod2vr) ip u (x, t) = ir (mod2vr) (13) 

respectively if v 2 > 1 and v 2 < 1. The former (resp. latter) is clearly 
stable (resp. unstable) because it corresponds to all pendula hanging 
downwards (resp. standing upwards) in the model of fig. ^ The same 
constant solutions arise from considering the constant g = {2k + 
corresponding to e = 62 = 1. 

2. If — 1 < e < 1 (e = ei in Fig. |2J) then the corresponding solution ffoO^e) 
can be written in terms of elliptic functions and describes the motion of 
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a particle confined in an interval contained in ] — tt, tt[ and oscillating 
around g = (mod27r) with some period H p (e). For <p, this translates 
into a periodic oscillating wave travelling with velocity v : If v 2 > 1 then 
tp oscillates around the stable equilibrium solution ip s = and describes 
a "plasma wave", see fig. 01(c), (d); if v 2 < 1 then tp oscillates around 
the unstable equilibrium solution tp w = it. Both kinds of tp are however 
unstable [El ED El- 

3. If e = &2 = 1, beside the constant solution yielding (|13[). there are in 
addition the solutions 

5o(±0 = 4tan- 1 [exp(±0]-^ 

Mod 2n, <7o(£) ~~ *■ i 71 " as £ — > ±00: the particle, confined in the interval 
—Tr<g<Tr, starts at 'time' £ = —00 from one top of the energy potential 
and reaches the other one at £ = 00. Mod. 2n, they translate into the 
following solutions of the original problem (jHJ) 



tp^(x,t;v) = 4 tan 1 |exp 
<Pq{x, t; v) = 4 tan" 1 |exp 



x — 
, x — vt 



vr 



7T if t> 2 > 1, 

if < 1. (14) 



Both families represent solutions (rotating clockwise or anti-clockwise ac- 
cording to the sign ±) with localized region of variation of tp and travelling 
with velocity v. The first families are clearly unstable because they de- 
scribe solutions of the model of fig. H where all pendula stand upwards 
except in the small moving region where they twist once around the axis; 
their values of the total mechanical energy H are infinite. The second 
families are stable |18l lTll3]. as it can be expected from the fact that they 
describe solutions of the mechanical model of fig. H where all pendula 
hang downwards except in the small moving region where they twist once 
around the axis, see fig. El (a); their values of the total mechanical energy 
H are finite. They represent respectively a soliton (tp^ ) and an antisoliton 
(<Pq) travelling with velocity v. Note that they fulfill © with n = ±1. 
For v = we have a static (anti)soliton. 

4. If e = e3 > 1 (see Fig. ^ then the corresponding solution go(±£;e) 
describes a particle moving towards the right and the left respectively, 
'for ever', since it has a sufficient energy to overcome the tops of the 
energy potential. Moreover, its kinetic energy and velocity are periodic 
with some period Ho(e), (with So — > 00 as e [ 1). This means, for any 
£ G R, 

3o(£ + H ;e) =5 (e;e) + 2vr, (15) 

i.e. is the sum of a linear and of a periodic function. Again, the 

corresponding solutions of the original problem © are unstable or ^3 ^ 
E] stable according to v 2 > 1 or v 2 < 1, because they correspond to 'most' 
pendula up or down in the model of fig. E The stable solutions (v 2 < 1) 
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([>Q describe evenly spaced "arrays of solitons and antisolitons" , travelling 
with velocity v, see fig. 01(b). 

No n with \n\ > 1 is allowed for travelling wave solutions fulfilling @. 

3 Adding the forcing and the first order dissipative 
terms in the equation 

We still adopt the Ansatz © for the solutions of and ask how the classes 
of solutions found in the preceding section are deformed. We are interested 
in stable solutions ip with bounded derivatives for x — > ±00. Having in mind 
the mechanical analog of the chain of pendula, it is natural to expect that 
solitonic solutions will survive perturbation: if we just smoothly switch 7 then 
even the static soliton will start to move and gradually accelerate; if we now 
also switch a, we expect that its propagation will approach a steady regime in 
which dissipation and forcing balance each other. We are actually going to see 
that not only the solitonic, but also the other classes of solutions found in the 
previous section survive perturbation whenever a compensation of the forcing 
term 7 with the dissipation term in a can take place. 
Replacing the Ansatz (j3J) in Q we find the equation 

(v 2 — l)g" — avg' + sing + 7 = 0, x G R. (16) 

If v = ±1 the second order derivative in Ijlfijl disappears and we get 

±ag = (sing + 7). 

Unless g is constant (and therefore equal to — sin -1 7 or 7r+sin _1 7) then there 
exists a £0 such that g'(£o) 7^ and go := g(Co) ~ sin _1 7, it + sin _1 7. 
Integrating in a neighbourhood of £0 one finds 

£ 9 

dz 



i-io = J di = ±a I 



sin z + 7 



Co 9o 

As g approaches respectively — sin -1 7 or -zr+sin -1 7 (mod. 2tt) the denomina- 
tor goes to zero linearly while keeping the same sign, and therefore the integral 
diverges logarithmically, implying that the corresponding time £ goes respec- 
tively to ±00 (or viceversa). The corresponding solution (p for Q is unstable, 
therefore is not interesting for our scopes, because: it yields a maximum of the 
energy density h as either £ — * 00, or £ — > —00, i.e. there it corresponds to 
all pendula standing upwards (while hanging downwards resp. as £ — > — 00 or 

i -> 00). 

If v 2 7^ 1 we perform the redefinitions [compare with 

e := -sign(^)^ 9(0--=-g(ti) if^ 2 >l, 

e-sign^)^ 0(0 Xv 2 <l, (17) 

£:=i=x 5(0 :=<?(£) -7T if« = 0; 
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Figure 3: Photographs of the mechanical model of fig. ^ (a) single soliton 
solution, (b) evenly spaced array of solitons, large (c) and small (d) amplitude 
plasma wave solutions (After A. C. Scott ^2], courtesy of A. Barone, see [2]) 
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Figure 4: The potential energy U for 7 = .1 and three mechanical-energy levels 
for fi = (no dissipation). 



we obtain 

g" + fig' + U g (g) = 0, (Gi, (18) 

which can be regarded as the 1-dimensional equation of motion w.r.t. the 'time' 
£ of a particle with unit mass, a 'potential energy' U{g) = —(cos g + 75) (see 
Fig. |1J) and a viscous force with viscosity coefficient 

a 

" ^ 7F^ (19) 

[in other words, in the equation a, v appear only through their combination 
1)191) ]. One immediately finds that the 'mechanical energy' is not 'time' inde- 
pendent, but decreases according to 



e' = -fig' 2 < 
U admits local minima (resp. maxima) in the points 

g™ := sin" 1 7+2fc7r (resp. gjf := tt — sin -1 7+2A;7r) 

and their corresponding values of U are 

W) = 



(20) 



7(sin 1 7 + 2kir) + \J\ — j 2 
7 (- sin" 1 7 + (2fc+l)7r) - \f\ 



r 



they linearly decrease with fc. 

We have to look for solutions of Q18|) such that 5' is bounded all over 
Multiplying Q18|) by 5' and integrating over an interval [£, £] we find 



AO = + ^ U(0 U(0 +m /Ae'K'- 



(21) 
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1. The constant (in £) solutions corresponding to the local minima, maxima 
of the potential energy become g(£) = g™, #(£) = gjf . The latter are the 
solutions for which there exists a £ such that <?(£) = and </(£) = 
(actually this happens for all £ € R). The particular ones <?(£) = p™, 
g(£) = gff correspond to the energy levels e = eo := U(g™), e = &2 '■= 
U{g^f) in Fig. ^ These solutions translate into the following constant 
solutions for the original problem: (mod. 2tt) 

ip s (x,t) = -shr 1 -/ ip u (x,t)= sin" 1 7 + vr (22) 

(the former are stable, the latter unstable). 

2. If there exists some 'time' £ such that g(£) = gjf and </(£) < (particle 
located at some maximum of {/ and moving leftwards), then necessarily 3 

#(£) — ► 00, 5 (£) — > -00, 

so for our scopes we can readily exclude this case. 

3. If there exists no 'time' £ such that g(£) is a local maximum of U, then 
necessarily the motion is confined in some interval ]s|!li><?j(f [ and |</(£)| 
by <|2~T|) keeps bounded. We have already treated the constant #(£) = gj™ 
solutions, so let us consider the others. 

If = (note that this is the case not only if a = 0, but also if a > 
and v = 0, i.e. for a static </?) and: 

(a) e < U(g^ 1 ), then <?(£) will be periodic, oscillating 'for ever' around 
<7™; in fig. 0]this corresponds to e=ei with g €]g-i, gff [■ The corre- 
sponding solutions ip of the original problem will describe again un- 
stable oscillations around the unstable equilibrium position if v 2 < 1, 
and unstable 'plasma wave' oscillations around the stable equilibrium 
position if v 2 > 1. 

(b) e = U(g^) , then #(£) -> c?f for both £ -> ±00; in fig. H this 
corresponds to e=e2 with g ^\g^x^9^\- If f 2 > 1 this will yield 
again an unstable solution tp, because the latter corresponds to all 
the pendula standing upwards at infinity, whereas if v 2 < 1 this will 
yield a new type of solution ip of the original problem, a kind of 
rigidly bounded soliton-antisoliton pair. Whether this is stable or 
not should be investigated. 

If /i (and therefore a) is positive, then necessarily 4 

m^=? 0, </(£) ^ 0. 

3 In fact, since U(g^) > U(g) for all g > gj^ , each term at the rhs of 12 II is positive for all 
£ < £, implying that |</(£)l > > for all £ < £. Hence the first limit in the previous 

relation. Now for g(£) growing without bound as £ — > —00 the local maxima of U^g(^)J 

decrease linearly, implying at least a linear growth of the rhs of 1211 1. whence the second limit. 
4 The second limit follows from 121 Hi and the existence of the lowest bound U{g™) for c, 
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As a consequence, using the equation of motion (|18|). as £ — > ±00 5 
must go to one of the following values: g^i, 9™ , 9k 1 ■ If in addition there 
exists some 'time' £ such that </(£) < (particle moving leftwards), then 
necessarily 5 

0(0 — ► 9k or 5fc-i, 5(0 — > 9k ■ 

Such a g translates into an unstable solution ip of the original problem, 
because it yields a maximum of the energy density h, i.e. it corresponds 
to all pendula standing upwards, for either £ — ► 00 or £ — ► —00 [depending 
which of the two possible definitions of g in 1)11(1 is adopted]. 

We are left with the last, most interesting possibility, namely a solution 
g such that 

5(0 — ► 9 k -i 9(0 — ► 9 k ( 23 ) 

and <?'(£) > for all £. This will translate into an unstable solution ip 
(almost all pendula standing upwards) if v 2 > 1 [see pi))], candidate 
stable solutions ip (almost all pendula hanging downwards) if v 2 < 1. 
The latter will describe the propagation of an (anti)soliton with velocity 
v [fig. |21 (a)]. That such a solution exists can be argued as follows. 
Impose just 1)23)1 1 . Consider first the case \i = 0: the mechanical energy 
of the particle has the constant value e = U(gjf_ l ), the corresponding 
solution will reach and overcome gjf after a finite time. Second, it is 
easily expected and not difficult to show |S] that for sufficiently large \x the 
corresponding solution will never reach gjf , but rather invert its motion 
at some g < gjf and go to g™ as £ — > 00. By continuity, we expect that 
there exists a special value p,(j) such that the corresponding solution g 
fulfills also (|23|) 9. As an immediate consequence of the definition of fi the 
propagation velocity v (for v 2 < 1) is determined [see formula 1)27)1 below]. 
One can determine ^((7) perturbatively in 7 at lowest order one finds 
p,("f) = 7T7/4 + which gives a v coinciding with (J2J), but higher orders 
will correct the latter formula. 

4. It remains to consider the cases that there exists some 'time' £& such that 
g(£k) = 9k (particle located at some maximum of U), but </(£&) > 

{ 

using standard methods in stability theory. As for the first limit, note that /(£) := J g' 2 (£,')d£,' 

( 

is nonnegative and monotonic. If per absurdum did not converge to a finite Im > but 
diverged as £ — > —00, then so would do the rhs l|21|l (because the term in square bracket is 
bounded for g G]ffk_i, the lhs and hence \g'\, in contradiction with the motion being 

confined in the interval. ^ — ?° Im implies a fortiori the first limit, as claimed. 

The first limit follows from excluding g™, which would imply the constant = g™ 
solution. As for the second limit, we can readily exclude gj^_i by the decrease law 1201 : 

we have to exclude g^ 1 by the same law, if g(£) 4 — ?° g^f , and by the existence of £, if 
9(0 — ► 9k-i- 
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(particle moving rightwards) at all such 'times'. Then necessarily </(£) > 
for all (finite) £ < As a consequence, it will be 



either <?(£) ^ — ?° — oo and ^ — 3 d oo 
or g(0 ^ M and f:= ^ oo 



(24) 



with some I < k, with |</(£)| remaining bounded. 

Consider first the case (|24[h . Employing as before a continuity argument 
in fj, and the invariance of Eq. (|18|) under o — > 17 + 27r, we actually infer 
that for any <?q > there exists a special value fi(^,g' ) and a finite time 
H(7,<7q) such that the corresponding solution <?(£;<?q,7) fulfills 

3(0; 3o, 7) =3^i 5'(0;3o,7)=3o 
g{Z; S / Q , 1 )=g*f 5'(S;3o,7)=3o 

and consequently, more generally, 

g(t + S) = 5(0 + 2vr. (25) 

Since, as one can expect and as we shall prove in [S], the dependence 
of S on (/g is strictly monotonic, one can choose H instead of g' Q as an 
independent variable. Again, this may translate into stable solutions Cp 
for Q only if < 1, and, as an immediate consequence of the definition 
of n, the propagation velocity v is determined [see formula (|2^)) below]. 
Therefore ip will describe an array of (anti)solitons travelling with such a 
velocity [fig. 01(d)]. The (anti)soliton solutions can also be regarded 
and obtained as g' j limits of (p. 

As for the cases (|24Jl o. we have hints [£j that the corresponding g will 
approach a g as £ —* 00. 

We collect our main results in 

Proposition 1 Mod. 2-7r ; stable, travelling-wave solutions of £7J) (where 1 > 
7 > and a > having bounded derivatives at infinity can be only of the 
following types (with £ := (x— vt)/y/\l — v 2 \): 

1. The static solution ip s (x,t) = — sin _1 7. 

2. The soliton/antisoliton solutions 0^(x,t) = g(±£) 

lim (x, t) = — sin -1 7 lim <p^(x, t) = — sin -1 7 ± 2ir (26) 



6 Otherwise, per absurdurn denote by | the largest £ < such that </(£) = 0. Then by 
1181 it is necessarily — U g I ) = g"(0 > 0. Since it cannot be U g I <?(£) I = 0, it must 



necessarily be g"(£) > 0, what implies </(£) < in a left neighbourhood of £ (namely £ is a 
time of inversion of the motion); by 1201 . going further backwards in time e will not decrease, 
and therefore there will be a £' k < £ < such that = gif and </(£*,) < 0. 
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travelling rightwards/leftwards resp. with velocity 



v = ± - A(7) ee ±v; (27) 
yfa 2 +/i 2 (7) 

at lowest order in 7 the function fi is given by ft{pf) = irj/4 + .... 

3. For any period S g]0,oo[ the "arrays of solitons/ antisolitons" solutions 
ip^ix^t) = g(±£), where 

g(Z + Z)=g(0 + 2Tr (28) 
travelling rightwards/leftwards resp. with velocity 

v = ±-^= = ±«. (29) 

jj, (as well as g itself) is a function 0/7,3. 
4- The "half-array of solitons/antisolitons" solutions <p > (x,i) = <7(±£) 

lim g(£) = — sin -1 7 lim g(£) = 00. (30) 

trevelling rightwards/leftwards resp. with velocity ±v. 

Remark 1. We emphasize that, in contrast with the unperturbed soliton 
(and array of solitons) solutions, where v was a free parameter (of modulus 
less than 1), v for the corresponding perturbed soliton (and array of solitons) 
solutions is predicted as the function of 7, a, g' given by formulae (|27|). (|29|). 

Remark 2. Note that g is a solution of 1)18(1 also if we identify £ + S with 
£, i.e. define £ (and therefore also x) as a point on a circle of length a multiple 
m5 of H, m = 1, 2, .... The corresponding solution ip of (JJ) will be defined for 
cc belonging to a circle of length mS\/l — u 2 , as well! 

Remark 3. The "half-array of (anti) solitons" solutions have no analog in 
the pure sine-Gordon case. 

Remark 4. In the list we may have to add for either a = and v 2 < 1, or 
a > and u = 0, the mentioned solutions tp(x, t) = <?(±£) such that 

lim (fix, t) = — sin" 1 7 = lim tp(x,t), (31) 

in case investigation should establish their stability. They would describe rigidly 
bounded soliton-antisoliton pairs moving with velocity v. 
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